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Abstract 

Several recent imaging experiments access the equilibrium density profiles of interacting particles 
confined to a two-dimensional substrate. When these particles are in a fluid phase, we show 
that such data yields precise information regarding substrate disorder as reflected in one-point 
functions and two-point correlations of the fluid. Using Monte Carlo simulations and replica 
generalizations of liquid state theories, we extract unusual two-point correlations of time-averaged 
density inhomogeneities induced by disorder. Distribution functions such as these have not hitherto 
been measured but should be experimentally accessible. 
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Place a perfect crystal in a disordered background deriving from a large number of ran- 
domly placed, quenched point pinning sites. Allow the crystal to relax to its minimum free 
energy state in this background. A qualitative picture of structure in the state which results 
is obtained by balancing the energy cost of distortions in the displacement field with the 
energy gained from accomodating locally to the pinning, a competition conventionally de- 
scribed via the paradigm of an elastic manifold in a random medium GO . Such a starting 
point, however, is clearly inadequate when equilibrium correlations in the pure system do 
not resemble those of a crystal or indeed of any system in which particles of fixed connectiv- 
ity are separated by a mean distance from which deviations are energetically penalized. The 
appropriate problem in this case is that of determining the correlations of an equilibrium 
fluid in a quenched random environment [3|] . Such a description applies to a large number of 
experimental systems. Indeed, systems modelled as elastic manifolds at low temperatures 
generically melt into fluid states at high enough temperatures. 

Equilibrium fluids in a quenched disordered background exhibit disorder-averaged ver- 
sions of the conventional correlation functions of the pure system. In addition, new classes of 
non-trivial correlations emerge 3|. Modern imaging techniques enable direct microscopic vi- 
sualization of several two-dimensional situations which may be mapped onto the problem of 
interest here: these include colloidal particles in their fluid phase on rough substrates 0,0,0], 
magnetic bubble domain arrays at large effective temperatures charge-density wave 
systems^] and varied phases of vortices in thin superconducting films P. Ilfl llll . I12I 
As we show in this Letter, such experiments raise the possibility of measuring novel cor- 
relation functions associated with fluid systems in the presence of quenched disorder. In 
addition, the topography of surface randomness can be accessed from the images yielded by 
such methods, even though it is the fluid (adsorbate) particles which are being imaged and 
not the substrate directly. 

This Letter describes a procedure for reconstituting a disorder potential in microscopic 
detail from knowledge of the time averaged density p(r) of the adsorbate in the liquid state 
and its spatial correlations. To this end, we also present accurate benchmarks for replica- 
based liquid state theories for correlations in disordered systems, comparing the predictions 
of these theories with direct simulations. Our proposals are illustrated using Monte Carlo 
computer experiments on a model two-dimensional fluid in a disordered background. Fig- 
ure E exhibits our main result: Fig^ (a) shows the random substrate potential used in our 
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FIG. 1: (color-online) (a) Color plot of the disorder potential /3Vd(x, y) used in our Monte Carlo 
simulation. The plot is a single realization of a Gaussian random potential of strength a 2 = .1 and 
correlation length £ = .12o"o. (b) The reconstituted potential (see text) with the measured time 
averaged density p(r) (Fig. 2 (a)) as input. The correspondence between (a) and (b) is close but 
can be systematically improved by averaging p(r) for longer times. The a 2 and £ obtained from 
(b) agree with the input values in (a) to within a few percent. 

Monte Carlo calculations, while Fig. H](b) exhibits the reconstituted potential obtained from 

n 

the induced p(r) , using a replica-based liquid state theory [3J . As these figures illustrate, to 
an accuracy limited only by the computational effort required for a well-averaged p(r), our 
calculation reconstitutes the imposed disorder potential. 

Our model system comprises N = 780 particles confined to two dimensions and interact- 
ing via a inverse-twelfth power pair potential. These particles also interact with a quenched 
one-body disorder potential Vd(r), modelled as a zero-mean Gaussian random field with spec- 
ified variance and short-ranged spatial correlations. The interaction part of the Hamiltonian 
is thus: Him = ^J2i<j(ff-) 12 + J2i v d{ri). We set e = I and er = 1, thus setting energy and 
length scales. We use periodic boundary conditions, benchmarking our calculations in the 
pure limit against earlier work [3]. 

A one-body random Gaussian disorder potential with zero mean value and exponentially 
decaying correlations, V<j(r), is constructed following a method proposed by Chudnovsky 
and Dickmanjlfil]. The variance (a 2 ) of the Gaussian distribution and its spatial correlation 
length (£) specify the potential. We introduce disorder with a suitably small spatial correla- 
tion (£ = 0.12) in units of the interparticle spacing. We study the system in the fluid regime 
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FIG. 2: (color-online) (a) Color plot of the configuration averaged density p(r) for a single disorder 
configuration (see Fig 1(a)) from a portion of our simulation cell. Plot (a) shows p(r) obtained 
from our Monte Carlo simulations using N = 780 particles. The averaging is over 1.41 x 10 5 
configurations each separated by 100 MCS. The average density po = 0.9. Plot (b) show the same 
quantity obtained from our theory for the corresponding disorder potential and po. 

for po = 0.05, 0.3 and 0.9, both with relatively weak (a 2 = 0.01) and with stronger (upto 
a 2 = 1) disorder. These disorder strengths are relatively moderate in comparison to the 
strength of the pair interaction, justifying perturbation theory in the disorder potential. We 
disorder average over five and ten disorder realisations for the weak and stronger disorder 
cases respectively, rejecting approximately the first 10 4 — 10 5 Monte Carlo Steps (MCS). We 
average over about 10 3 — 10 5 configurations each separated by 10 2 MCS, ensuring adequate 
thermal averaging of the local density and correlation functions. 

For an equilibrium fluid in the absence of disorder, all points in space are equivalent. 
As a consequence, one-body distribution functions are structureless. This property does 
not hold, of course, within a particular realization of disorder. It is only restored upon a 
disorder average. If the fluid particles do not interact, the Boltzmann relation connects the 
time averaged density (p(r)) to the local potential via (p(r)) ~ p exp[— (3Vd(r)] where p 
is the average density of the fluid, (3 = l/k B T with T the temperature, and Vd(r) is the 
disorder potential. This result is valid only in the limit of vanishingly small po or as T — > oo; 
equivalently, in the zero correlation limit. Increasing p enhances the magnitude of the 
signal, making it easier to observe, but at the same time introduces non-trivial correlations 
in particle positions. Thus, in a given disorder configuration, one-particle distributions 
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reflect both the inhomogeneities of the potential as well as the consequent structuring of the 
local density field as a consequence of correlations [l^. 

Figs. 2 (a) and (b) show the time averaged density p(r) averaged over 1.41 x 10 5 indepen- 
dent configurations for p = .9 and disorder strength cr 2 = 0.1. Note the strong structuring 
visible in Fig 2(a). There are pronounced, correlated peaks and troughs in the time-averaged 
local density, indicative both of pinning due to local potential minima as well as the indirect 
effects of inter-particle correlations. Figs 2(b) shows the theoretically obtained plot of p(r) 
using a method described below. As can be seen, the theory captures the essential features 
of the simulation data with fair accuracy. 

An inhomogeneous potential V^(r) couples to the local number density and can thus be 
absorbed into the definition of the chemical potential /i, leading to p — > p(r) = \i + Sp(r) 
where 5p(r) = — V^(r). The Ursell function connects the local time-averaged density (p(r)) 
to the inhomogeneity in the chemical potential via 

(p(r)> = (p (r)> + P J dr'S(r, r')fy(r') + (1) 



where S(r, r') is defined as S(r,r f ) = (p(r)p(r')) — (p(r))(p(r'))[17|. Results perturbative in 
weak disorder are obtained by expanding about the pure limit, in which case (po(r)) = po, 
the average fluid density. The density p(r) can be generated if the appropriate correlations 
(p(r)p(r')) and (p(r))(p(r')) are available. To lowest order these are the correlations of the 
pure system. The results can be extended to larger disorder (a non-trivial initial (po(r))) 
provided accurate values of S(r, r') computed in the disordered background are available. 

For weak disorder, we may approximate (p(r)p(r')) — (p(r)(p(r')) by PoM|r ~~ r 'l) where 
h(r) is the pair correlation function of the pure system. At somewhat stronger disorder, an 
alternative approach improves on this result by using the disorder renormalized version of 
these correlation functions, replacing h{r) above by the function g^'(r) — g^ 2 \r) + <5(r)/po. 
Here g^'(r), the disorder averaged analog of the radial distribution function, is 

g {1) (r) = [(p(0)p(r))}/pl-5(r)/p . (2) 

The analog of an Edwards- Anderson parameter reflecting the correlations of time-averaged 
density inhomogeneities is the "off-diagonal" distribution function g( 2 \r), defined through 

9 (2) (r) = [(p(0)}(p(r))}/pl. (3) 
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Here, (...) denotes a thermal average for the disordered system prior to the disorder aver- 
aging, while the brackets [. . . ] denote an average over disorder. 

Thus, in the presence of correlations, the response of the density to a nonzero external 



potential (apart from an overall normalization) is given by 17(, 

<p(r)> = po - Ppl J dr' (g^(\ r - r' |) - g^(\ r - r' |)) ^(r') - (3p V d (r) + ... (4) 

At this point apart from specifying how g^(r) and g^ 2 \r) are obtained which we do 
below, our procedure for obtaining p(r) is complete. Using Eq. 0]we may obtain the density 
(Fig. 2 (b)) from the potential (Fig. 1(b)) or, by using a Fourier transform to invert Eq. 
obtain the disorder potential (Fig. 1(b)) from the "experimental" density (Fig. 2(a)). In 
principle, inversion of the density to obtain V^(r) requires the specification of g^ 2 \r) which 
itself depends on V^(r). In practice, a simple iteration starting from pure system correlations 
converges rapidly. 

Our calculation of g^'(r) and g^ 2 \r) is based on an early replica-based approach to the 
calculation of correlation functions in a disordered fluid 3). The replica method is applied 
to the partition function of a system of classical particles interacting via the Hamiltonian 
H = H kinetic + \Y J i^j V{\ r i~ r j + Vd{ri), with V(r) is a two-body interaction potential 
between the particles and V^(r) a quenched, random, one-body potential drawn from a 
Gaussian distribution of zero mean and short ranged correlations: [Vd(r)Vd(r')] = K{\ r — r' | 
) . One then notes that the replicated partition function resembles the partition function of a 
system of n "species" of particles, each labeled by an appropriate replica index, interacting 
via a two-body interaction which depends both on particle coordinates (r.j, i\,) and replica 
indices (a, 13). This system of n species of particles is considered to be a n-component 
mixture [l7|. Taking the n — >• limit in the Ornzstein-Zernike equations governing the 
properties of the mixture and assuming replica symmetry, yields the following coupled set 
of equations, in Fourier space, coupling the pair correlation functions and h^> to the 
appropriate direct correlation functions and C^ 2 \ 

(1) _ CW(k)-[CM(k)-CV(k)? (5) 



[1 - CM(k) + C^\k)] 2 ' 
The Fourier transforms involved are defined as <p(k) = Po j dr<p(r) exp (— ih ■ r). The 
replicated Ornstein-Zernike relations must be supplemented with specific closures. Previous 



work applied to the problem of flux-lattice melting jl8{] in the presence of quenched point 
pinning used the simplest such closure, the Hyper- Netted Chain (HNC) closure^. We have 
experimented with a variety of closure schemes to test out the accuracy of this approach 
to the calculation of fluid correlations. Our best results are obtained with a closure scheme 
due to Rogers and Young (RY) [19| for pure systems: 

exp [Y (1) (r)/(r)] - 1" 



C (1) (r) = exp[-(3(V {1 \r) + V {2 \r))] 
C<$(r) = -(3V {2 \r), 



1 + 



f(r) 



Y (1) (r), (6) 



where the function f(r) = 1 — exp(— ar) interpolates between HNC and Percus-Yevick (PY) 
values for large (small) r and a the "switching" parameter is chosen so as to enforce ther- 
modynamic consistency. Here Y^ v \r) = h^ u '(r) — C^(r), v — 1,2. The pair interaction 
potential is V^'if) and the off-diagonal potential V^(r) is the disorder averaged correla- 
tion function of the disorder potential V^(r). The integral equations for the liquid state 
correlations as given above must be solved numerically 3l. Il8l l20j]. 

The correlation function g^\r) obtained from self-consistent solutions of the integral 
equations is qualitatively similar to that for the pure case at low and intermediate levels 
of disorder. At large values of disorder, a disorder- induced suppression of structure in 
gW(r) is clearly apparent. This is shown in Fig. 3(a) which exhibits the quantity Ag(r) = 
grW(r) — g(r), the difference between the pair distribution function in the disordered and 
pure cases. A direct comparison to the simulation data is shown in Fig. 3 (b). We find that 
while HNC and PY respectively under and over estimate Uj correlations, the RY closure 
gives very good agreement with our simulations. (This is consistent with results for the 
pure system.) The parameter a in the latter is taken to be its pure system value of 0.3, a 
value which agrees best with simulation data over the full range of densities and disorder 
strengths. 

We show the off-diagonal correlation function g^ 2 '(r) at a 2 = 0.01 for p — -05, .3 and 
.9 in Fig. 4(a). In the simulations, this correlation function is computed by extracting 
the correlations of the local time-averaged density in a particular configuration of disorder 
and then averaging over several disorder realizations. Correlations increase with increasing 
density as expected, as manifest in the oscillations of g( 2 \r). For the same reason, the off- 
diagonal correlations for any density also increase with the strength of the disorder correlator 
(Fig. 4(b)). In both cases, our theoretical estimates based on the RY closure agree to good 
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FIG. 3: (a) Plot of Ag^(r) for a single density po = 0.9 and a 2 = 1.0. (b) Comparison of g^(r) 
for density po = 0.9 and a 2 = 1.0 of simulation data (MC) with the results of HNC, PY and RY 
closures. Ten disorder configurations were used in the averaging. 
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FIG. 4: (a) Plot of g^(r) comparing simulations (points) and result from the RY closure (lines) 
for densities po = -05 (+),.3(o) and .9(box) for a 2 = .01. Ten disorder configurations were used 
in the averaging, (b) Plot of g^ 2 \r) - simulations and theory (RY) for a 2 = .01 (+), .1 (o) and 1 
(box) for po = 0.9 

accuracy with correlation functions obtained from the simulations. 

Liquid state theory based approaches are easily adapted to a variety of different problems 
vis a vis. direct simulations. Thus, the approach presented here is potentially useful to a 
variety of theoretical models for experimental data. The following limitations must, however, 
be kept in mind: First, p(r) must be accurately determined from experimental pictures for 
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correlations to be obtained correctly. Since the data are often thresholded, reflecting some 
level of coarse graining of the amplitude of the signal, a direct comparison with theory 
may not always be straighforward. (Interestingly, the size of the probe particles do not 



impose any theoretical limit 



21 1 on the spatial resolution of V^(r) - this is limited only by 



the accuracy with which p(r) can be imaged.) Second, our analysis neglects higher order 
density-density correlations. Finally, at high disorder strengths and strong correlations, the 
possibility exists of a transition into a state in which replica symmetry is broken, a transition 



22]. 



which would not be captured in this calculation 

What types of experiments might access the unusual "off-diagonal" correlation function 
? Any experiment in which repeated 'snapshots" of the system are taken, with suffi- 
cient statistics, should yield results which can be addressed by these methods. For vortices 
in superconductors, scanning tunneling spectroscopy 



Lorentz and magnetic force microscopy 
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23 



24j magneto-optic imaging 11( and 



261 ] can all be used as structural probes of vortex 



structure in thin superconducting films. Magnetic bubble domain arrays are another exper- 
imental system in which many of the features described here should be accessible fZj. We 
also note that several recent experiments directly access the configurations of a large number 
of colloidal particles confined to two dimensions I28. 29, 30]. In these experiments, ex- 
tensive statistics for particle positions can be generated. The colloid literature has hitherto 
concentrated on making substrates as homogeneous and disorder-free as possible. However, 
as we argue here, interacting colloidal particles moving on disordered substrates can ex- 
hibit interesting and non-trivial correlations which have not hitherto been characterized. 
Experimental work in this direction would be welcome. 
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